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Abstract 

We study the earth matter effects in the Tokai-to-Kamioka-and-Korea exper- 
iment (T2KK), which is a proposed extension of the T2K (Tokai-to-Kamioka) 
neutrino oscillation experiment between J-PARC at Tokai and Super-Kamiokande 
(SK) in Kamioka, where an additional detector is placed in Korea along the same 
neutrino beam line. By using recent geophysical measurements, we examine the 
earth matter effects on the oscillation probabilities at Kamioka and Korea. The 
r-^. ' average matter density along the Tokai-to-Kamioka baseline is found to be 2.6 

g/cm , and that for the Tokai-to-Korea baseline is 2.85, 2.98, and 3.05 g/cm for 
the baseline length of L- = 1000, 1100, and 1200 km, respectively. The uncertainty 
of the average density is about 6%, which is determined by the uncertainty in the 
correlation between the accurately measured sound velocity and the matter den- 
sity. The effect of the matter density distribution along the baseline is studied by 
using the step function approximation and the Fourier analysis. We find that the 
V[_i — > v e oscillation probability is dictated mainly by the average matter density, 
with small but non-negligible contribution from the real part of the first Fourier 
mode. We also find that the sensitivity of the T2KK experiment on the neutrino 
mass hierarchy does not improve significantly by reducing the matter density error 
from 6% to 3%, since the measurement is limited by statistics for the minimum 
scenario of T2KK with SK at Kamioka and a 100 kt detector in Korea considered 
in this report. The sensitivity of the T2KK experiment on the neutrino mass hi- 
erarchy improves significantly by splitting the total beam time into neutrino and 
anti-neutrino runs, because the matter effect term contributes to the oscillation 
amplitudes with the opposite sign. 
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1 Introduction 

Many neutrino oscillation experiments have been performed to measure the neutrino 
oscillation parameters. In the three neutrino model, there are nine fundamental param- 
eters; three masses (mi, m 2 , 7773), three mixing angles (#12, #13, #23) and three phases 
(5, (pi, (f>2) of the lepton-flavor- mixing matrix, MNS matrix [T|. Among the three CP 
phases, one (5) is lepton number conserving and two (0i, <p2) are lepton number non- 
conserving phases of the neutrino Majorana masses. Out of the nine parameters, the 
neutrino-flavor oscillation experiments can measure six parameters; two mass-squared 
differences {m\ — m\, m\ — m\), all the three mixing angles (#12, #13, #23), and the lepton 
number conserving phase (5). Among them, the magnitude of the larger mass-squared 
difference, \m\ — m\\, and a mixing angle, #23 , have been measured by the atmospheric 
neutrino experiments [21 [3], which have been confirmed by the accelerator based long- 
baseline (LBL) neutrino oscillation experiments, K2K [I] and MINOS [5]. The smaller 
mass-squared difference, m\ — m 2 , and the corresponding mixing angle, #12, have been 
measured by the solar neutrino oscillation experiments [6], and the KamLAND exper- 
iment [7] that observes the oscillation of reactor anti-neutrinos at more than 100 km 
distance. The third mixing angle, #13, has been looked for in the oscillation of reactor 
anti-neutrinos at about 1km distance. No oscillation has been observed and only the 
upper bound on sin 2 20i 3 of about 0.14 has been reported [8] [9]. 

Summing up, out of the nine parameters of the three neutrino model, six parameters 
can be observed in neutrino-flavor oscillation experiments, among which \m\ — m 2 |, #23, 
m 2 — m i' an d $12 have been measured, and the upper bound on sin 2 26 , 13 has been 
obtained. The goal of the next generation neutrino oscillation experiments is hence 
to determine the remaining parameters, i.e., the sign of rnj — m 2 (the neutrino mass 
hierarchy pattern), a finite value of the third mixing angle #13, and the CP violating 
phase of the lepton sector, 5. 

Several neutrino oscillation experiments are starting or being prepared aiming at the 
measurements of #13 , which can be categorized into two types; one is the next generation 
of reactor neutrino experiments, and the others are accelerator based neutrino oscillation 
experiments. Three experiments, Double CHOOZ [10], RENO [UJ, and Daya Bay|12j. 
will be able to observe the oscillation (disappearance) of z/ e 's from reactors at about 1km 
distances if sin 2 2# 13 £ 0.01. In accelerator based LBL experiments, v^ or v^ beam from 
high energy tt + or 7r~ decay- in-flight, respectively, will be detected at a distance of a few 
to several hundred km's away and the v^ — > v e or z/ M — > v e transition is looked for. The 
transition rate is proportional to sin 2 #13, and the next generation LBL experiments, the 
Tokai-to-Kamioka (T2K) experiment [13] and the NO^A experiment [H], have a chance 
to discover the v^ — > v e transition, if sin 2 26*13 ~ 0.02 1 . 

In this report, we focus our attention on the proposed one-beam two-detector ex- 
periment, the Tokai-to-Kamioka-and-Korea (T2KK) experiment [T5l |16| \TT \ [T8 | [i~9 j 120]. 
T2KK is an extension of the T2K experiment where an additional huge detector is 
placed in Korea along the T2K neutrino beam line between J-PARC at Tokai and Super- 
Kamiokande (SK) at Kamioka. We show the cross section view of the T2KK experiment 
in Fig. [TJ The center of the T2K neutrino beam has been designed to go through the 
underground beneath SK so that an off-axis beam (OAB) between 2° and 3° upward 



1 Recently, the T2K experiment [37] and the MINOS experiment [55] reported hints of non-zero #is 
based on observation of v u — > v e candidate events. 
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Figure 1: Schematic cross section view of the T2KK experiment when the center of the 
neutrino beam from J-PARC is 3° below the Super Kamiokande detector. The center 
of the beam is shown by the red arrows, and the black arrows show the baselines of a 
few off-axis beams (OAB). 



from the beam center can be observed at SK, in order to maximize the sensitivity to the 
V/j, — > v e oscillation. It then follows that the center of the T2K beam reaches the sea level 
in the Japan/East sea, and the lower-side of the same beam at 0.5° or larger off-axis 
angle goes through Korea [151 E] • Therefore, if we place an additional neutrino detector 
in Korea along the T2K beam line, we can perform two LBL experiments at the same 
time. It has been demonstrated that the T2KK experiment has an outstanding physics 
potential to determine the neutrino mass hierarchy pattern, normal (m\ — mf > 0) or 
inverted {m\ — m\ < 0), and to determine the CP violating phase 5 [16j [T7J \TE\ . 

The key of determining the neutrino mass hierarchy pattern in the T2KK experiment 
is the earth matter effect, which arises from the coherent interaction of v e (or z/ e ) off 
the electrons in the earth matter [2T]. The matter effect enhances or suppresses the 
Vp, —?• v e oscillation probability, respectively, if the mass hierarchy pattern is normal 
(mf — m 2 > 0) or inverted (m| — m\ < 0). Furthermore, the matter effect around the 
first oscillation maximum in Korea (L £ 1000 km) is significantly larger than the one at 
SK (L = 295 km) because its magnitude grows with neutrino energy, and the oscillation 
phase is roughly proportional to the ratio L/E. Fortunately, we can observe the first 
oscillation maximum at both detectors efficiently, if a detector is placed in the east coast 
of Korea such that 0.5° ~ 1° OAB can be observed at L ~ 1000 km when the beam 
is adjusted to give the 3.0° ~ 2.5° OAB at SK pU [18]. In a previous work [H], it 
has been found that the neutrino mass hierarchy can be determined most effectively by 
observing 0.5° off-axis beam in Korea at L ~ 1000 km and 3.0° off-axis beam at Kamioka. 
With this combination, the neutrino mass hierarchy pattern can be determined at 3-cr 
level when sin 2 2#i3 £ 0.06, if we place a 100 kt level Water Cerenkov detector in Korea 
during the T2K experimental period with 5 x 10 21 POT (protons on target). This is the 
minimum scenario where no additional detector is placed in Kamioka and no increase is 
the J-PARC beam power is assumed. A 100 kt level fiducial volume for the detector in 
Korea is necessary in order to balance the statistics of SK with 20 kt at 1/3 the distance 
from the J-PARC. A more grandiose scenario with huge detector both at Kamioka and 
Korea has been considered in Ref. |16j . 

Since the earth matter effect is very important in estimating the physics discovery 
potential of the T2KK proposal, we study in this report its impacts and uncertainty in 
detail by using the recent geophysical measurements of the earth matter density beneath 



Japan, Japan/East sea, and Korea. Both the average and the distribution of the matter 
density along the Tokai-to-Kamioka and the Tokai-to-Korea baselines are studied, as 
well as their errors. We also examine the impact of using both of the neutrino and anti- 
neutrino beams, since the earth matter effect on the v^ — > v e and v^ — > v e transitions 
have opposite signs. 

This paper is organized as follows. In section 2, we study the earth matter density 
profile along the Tokai-to-Kamioka and Tokai-to-Korea baselines, and then we evaluate 
the average and the Fourier modes of the matter densities along the baselines. In section 
3, we discuss how to treat neutrino oscillations with non- uniform matter distributions. 
We obtain both the exact solution to the oscillation amplitudes and on approximate 
analytic expression for the oscillation probability that makes use of the Fourier series 
of the matter distribution along the baseline. In section 4, we introduce our analysis 
method that quantifies the sensitivity of the T2KK experiment on the neutrino mass 
hierarchy pattern by introducing a x 2 function that accounts for both statistical and 
systematic errors. In section 5, we discuss the matter effect dependence of the capability 
of determining the mass hierarchy pattern in the T2KK experiment. In section 6, we 
consider the matter effect for the anti-neutrino experiment, and examine the impact of 
combining both neutrino and anti-neutrino oscillation experiments. Finally, we sum- 
maries our findings in section 7. In Appendix, we give the second-order perturbation 
formula for the neutrino oscillation amplitudes in terms of the matter effect term and 
the smaller mass squared difference. 

2 Earth matter profile along the T2KK baselines 

In this section, we study the matter density distribution along the baselines of Tokai- 
to-Kamioka and Tokai-to-Korea, and then estimate the average and the Fourier modes 
of the matter profile along the baselines. 

2.1 Tokai-to-Kamioka baseline 

First, we show the Tokai-to-Kamioka (T2K) baseline and its cross section view along 
the baseline in Fig. |2j In the upper figure, the scale along the baseline shows the distance 
from J-PARC by tick marks every 50 km. The shaded region is called Fossa Magna [22J , 
which is the border between the two continental plates, North American Plate and 
Eurasian Plate. In this region, the sediment layer is estimated as 6 times deeper than 
that in the surrounding region. The lower plot shows the cross section view of the T2K 
experiment. The horizontal line shows the distance from the J-PARC and the vertical 
axis gives the depth from the sea level. The number in each region stands for the average 
matter density in units of g/cm . In this figure, the thickness of the sediment layer except 
for the Fossa Magna region is assumed to be 1km [23J and we refer to the geological map 
in Ref . [21] . In Fossa Magna, the sediment layer is as deep as 6 km [22] , which is mainly 
composed of limestone and sandstone with others [21| . According to Ref. [22] , the density 
of such sediment layer is 2.5g/cm . On the other hand, the sediment layer near Kamioka, 
L > 230 km along the baseline, is mainly composed of granite whose density is about 
2.8g/cm . 

When we average out the matter density along the baseline shown in Fig. 121 we obtain 
the average density of 2.6 g/cm , which is significantly lower than the value 2.8g/cm 
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Figure 2: Upper: The Tokai-to-Kamioka (T2K) baseline. The distance from J-PARC 
along the baseline is shown by tick marks every 50 km. The shaded region is called 
Fossa Magna [22], where the North American plate and the Eurasian plate meet. Lower: 
The cross section view of the T2K experiment. The horizontal line shows the distance 
from the J-PARC [25] and the vertical axis measures the depth from the sea level. The 
number in each region shows the mean matter density in units of g/cm . The neutrino 
beam line between J-PARC and Kamioka is shown by an arc. 
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Figure 3: The relation between the matter density and the p-wave sound velocity (V p ) 
according to Ref . [26] . The solid line gives the mean value of the density, and the dotted 
lines in both sides show the error of the estimated density. The dashed lines show how 
we estimate the density and its error when V p = 6.5 km/sec. The three square points 
are the reference points in Preliminarily Reference Earth Model (PREM) 



quoted in the Letter of Intent (LOI) of the T2K experiment |13j . The difference is mainly 
due to Fossa Magna, whose lower density has not been taken into account in the past. 

The error of the average density can be estimated from the error of the mean density 
in each region, and the uncertainty of the boundary of each layer. In most of research 
works in geophysics, these two values are evaluated by using the seismic wave observation. 
The mean value of the matter density is estimated from the velocity of the seismic wave 
in each region. In this work, we adopt the density-velocity correlation of Ref. [25] to 
estimate the matter density from the p-wave velocity. In Fig. [21 we show the density 
velocity relation. The solid line shows the relation between the sound velocity and the 
mean matter density according to Ref. [26J, which can be expressed as 



p = -0.00283V; 4 + 0.0704V"; - 0.598V"; + 2.23V P - 0.7 . 



(1) 



where V p is the p-wave sound velocity inside the matter in units of km/sec, and the 
matter density p is in units of g/cm . The dotted lines in both sides show the error 
of the estimated density, which is about 6%. The dashed lines show how we estimate 
the density and its error when V p = 6.5 km/sec. The three squares are the reference 
points shown in Preliminarily Reference Earth Model (PREM) [27], which is often used 
to estimate the average density along the baseline for various LBL experiments. These 
points lie close enough to the mean density line of eq. (0Q), confirming the consistency 
between the PREM and the other geological measurements. On the other hand, the 
location of the boundary is evaluated from the reflection point of the seismic wave. The 
error of the boundary depth between the sediment layer and the upper crust is estimated 
as ±300 m. It affects the average density along the T2K baseline only by 0.05%, and it 
can be safely neglected. The error of the assumed sound velocity measurement at various 
location can also be neglected. In this report, we assume that the matter density error 
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Figure 4: Fourier modes of the matter density distribution along the Tokai-to-Kamioka 
baseline in units of g/cm . Black circles and white squares show the real and the imagi- 
nary parts of the Fourier coefficients, respectively. 

is determined by the 6% systematic uncertainty in the density-velocity relation, which is 
taken to be common (100% correlation) at all location along the baseline. On the other 
hand, the error of the average matter density along the T2K baseline and the Tokai- 
to-Korea baselines are taken to be independent in order to make the most conservative 
estimate. The 6% error adapted in this report is a factor of two larger than the error 
adapted in Refs. [TT1 IT8] . The matter density error is expected to become smaller by 
making use of the information from the other measurements, e.g. the gravity anomaly, 
the magnetic anomaly, and by actually digging into the crust. 

We note that the density fluctuation along the T2K baseline of Fig.[2]from the average 
of 2.6 g/cm is about 8%, which is comparable to the 6% error of the average matter 
density It is therefore of our concern that the approximation of using only the average 
matter density along the baseline may not be accurate enough. In Ref. [23], the authors 
have shown that the Fourier analysis is useful to study the impacts of the matter density 
distribution. We therefore show the Fourier coefficients of the matter density profile 
along the Tokai-to-Kamioka baseline in Fig. HI Black circles and white squares show the 
real and the imaginary parts of the Fourier modes in units of g/cm , respectively, which 
are defined as 
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where p(x) is the matter density along the baseline at a distance x from J- PARC. Positive 
value of Re(pi) is clearly due to the low density of Fossa Magna in the middle of the 



baseline, whereas the negative values of Im(/?2) and Im(/>3) reflect the high matter density 
in the Kamioka area that makes the distribution asymmetric about the half point (x = 
L/2) of the baseline. We note, however, that the magnitude of all the Fourier modes is less 
than 0.05 g/cm , which is less than 2% of the average matter density of po = 2.6g/cm . 
We can therefore expect that the fluctuation effect can be safely neglected in the T2K 
experiment, as will be verified below. 

2.2 Tokai-to-Korea baselines 

Let us now examine the matter density distribution along the Tokai-to-Korea baseline. 
In Fig. [5j we show the surface map of the T2KK experiment (up), its cross section view 
(middle), and the matter density distribution along the baselines (bottom). Along the 
center of the neutrino beam line shown in the surface map, the distance from J-PARC 
is given by tick marks every 100 km. The areas enclosed by the three circles show the 
regions which have been studied by geophysicists; Yamato basin, Oki trough [28J and 
Tsushima/Ulleung basin [2HI ED]- In the middle figure, the five curves correspond to the 
baselines for L = 1000, 1050, 1100, 1150 and 1200 km, where the horizontal axis measures 
the distance along the earth arc from J-PARC, and the vertical axis measures the location 
of the baselines below the sea level. The numbers in white squares represent the average 
matter density, in units of g/cm 3 , of the surrounding stratum: the sediment layer (2.5), 
the upper crust (2.8), the lower crust (2.9), and the upper mantle (3.3). A lump of 
3.1 g/cm at Korean coast is a lava boulder [30J. In order to produce this figure, we 
refer to the recent study on the Conrad discontinuity, the boundary between the upper 
and the lower crust, and on the Mohorovicic discontinuity (Moho-discontinuity), the 
boundary between the crust and the mantle, below Japan and Korea [3T] [H2]- Regarding 
the matter profile beneath the Japan/East sea, we divide the sea into three regions, 
Oki trough (east side) [28], Oki island (middle) [29], and Tsushima/Ulleung basin (west 
side)[30j. According to Refs. [281 E0], typical depths of the Conrad and the Moho- 
discontinuities below Oki trough (Tsushima/Ulleung basin) are 8 km (7km) and 20 km 
(17 km), respectively. On the other hand, Ref. [29J shows that they are 6km and 25 
km, respectively, below the sea level at around Oki island. Although there is no direct 
measurement of the discontinuities along the Tokai-to-Korea baselines, we may assume 
that the depth of the discontinuities is roughly characterized by the depth of the seabed 
of each region. If the sea is shallower, the discontinuities should be deeper, and vice 
versa. Because the seabed above the baseline is rather flat the east and west sides of Oki 
islands, we obtain the estimates of Fig. E]from the measured depths of the discontinuities 
in each region. 

In addition, we adopt the following simple interpolations around the sea bounds. As 
for the Moho-discontinuity, we assume that it is connected from 32 km below the sea level 
at the Japanese coast, the point (380 km, 32 km) in the middle figure, to 20 km below 
the Japan/East sea at (440 km, 20 km) by a straight line as illustrated in the figure. It 
is at this point, about 60 km west from the coast along the beam direction, the depth 
of the Japan/East sea reaches 1 km, the average depth of Oki trough. As for the Oki 
island region, the Moho-discontinuity is connected by straight line from (600 km, 20 km) 
to (630 km, 25 km), flat from 630 km to 700 km, and then by another straight line from 
(700km, 25km) to (730km, 17 km). Near the Korean coast, it is connected by a straight 
line from (970km, 17km) to (990km, 32 km). Here the slope is shaper because the Korean 
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Figure 5: The surface map (up) of the Tokai-to-Kamioka-and-Korea (T2KK) experi- 
ment, the cross section view (middle) and the density distribution along the baselines 
(bottom). The center of the neutrino beam line is shown in the upper figure, along which 
the distance from J-PARC is given by the tick marks every 100 km. The three circles 
show the regions which have been studied by geophysicists; Yamato basin, Oki trough 
[28] and Tsushima/Ulleung basin [29j [30]. In the middle figure, the five curves show 
the baselines for L = 1000, 1050, 1100, 1150 and 1200 km, where the horizontal axis 
gives the distance from J-PARC along the arc of the earth surface, and the vertical scale 
measures the location of the baseline below the sea level. The numbers in white squares 
represent the average matter density of the surrounding stratum, in units of g/cm 3 . 



coast is only about 20 km away from the western edge of Tsushima/Ulleung basin [30J . 
As for the Conrad discontinuity, we also assume the linear interpolation between (380 
km, 15 km) and (440 km, 8 km) at the Japanese coast, (600km, 8km) and (630km, 5km) 
at east side of Oki island region, (700km, 5km) and (750km, 7km) at west side of Oki 
island region, and between (970 km, 7 km) and (990 km, 17 km) at the Korean coast. 
The lava bounder lies on top of the oblique section of the Moho-discontinuity and it 
is assumed to have the triangle shape along the cross section with the apexes at, (970 
km, 17 km), (990km, 15km), and (990km, 30km) as also illustrated in the middle figure. 
Although the above treatments of the boundary conditions are very rough, we do not 
elaborate them further because the neutrino oscillation probabilities are very insensitive 
to the exact location of the boundaries. 

With the above setting, we obtain the matter density distribution along the baselines 
shown in the bottom figure of Fig. [5j Here the horizontal axis measures the distance x 
from J-PARC along the baseline. We assume that a detector in Korea is placed at the 
sea level (zero altitude) along the plane of the center and the center of the earth, at the 
distance L from J-PARC. Along the baseline of length L, the point x (km) away from 
J-PARC appears on the cross section figure in the middle of Fig. [5] at the distance I (km) 
along the earth arc, and the depth d (km) below the sea level; 



, x I R 2 -L 2 /A 

I = R arcsin \ —\ 



R\ R? -xL + x 2 



d = R _ v / x 2_ xL + R 2 = x ( L - x ) h+o(R- 3 )] , (3b) 

2R l V /J 

where R is the radius of the earth, 6378.1 km. The five arcs are the cross section view 
of Fig. [5] are then obtained parametrically as (l(x),d(x)). It is also noted that I is not 
much different from x for L ~ 1000 km; I is about 0.01% shorter than x at ~ 300 km, 
and 0.1% longer x at x ~ 1000 km. The effect of the lava boulder near the Korean coast 
is seen in the matter density distribution for L = 1100 km and 1150 km. We also notice 
that the baseline for L = 1000 km goes through the sediment layer (2.5 g/cm ) near the 
Korean east at 970km £ I £ 990km (969km £ x £ 989km). All the baselines go through 
the sediment layer near the J-PARC until the baseline reaches the depth of d = 1 km 
(0 < x < 10 ~ 13 km). 

Let us now estimate the average matter density and its error along each baseline. 
Since the traveling distance in the mantle and that in the crust depend on the locations 
of their boundaries, the error of the average matter density is affected by the errors of 
the depths of the discontinuities. It is especially important for the Moho-discontinuity 
because the density difference between the lower crust and the upper mantle is as large 
as 0.4g/cm , which is about 13% of their average. Therefore, we examine the impacts of 
varying the depth of the Moho-discontinuity on the average matter density. In Table HJ 
we show the average matter density along the baselines of L = 1000, 1050, 1100, 1150 and 
1200 km, in units of g/cm , for five locations of the Moho-discontinuity. The mean corre- 
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Baseline length 


The depth of the Moho-discontinuity below the sea level 


L 


-3a 


-la 


mean 


+la 
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1000 km 


2.923 


2.900 


2.854 


2.854 


2.854 


1050 km 


2.973 


2.954 


2.944 


2.933 


2.862 


1100 km 


3.001 


2.988 


2.980 


2.971 


2.952 


1150 km 


3.049 


3.039 


3.027 


3.002 


2.986 


1200 km 


3.054 


3.052 


3.051 


3.049 


3.036 


Repi (L = 1000 km) 


-0.097 


-0.076 


-0.031 


-0.031 


-0.031 


Repi (L = 1050 km) 


-0.102 


-0.098 


-0.094 


-0.090 


-0.028 


Repi (L = 1100 km) 


-0.117 


-0.110 


-0.101 


-0.105 


-0.099 



Table 1: Average matter density along the baselines of L — 1000, 1050, 1100, 1150 and 
1200 km, in units of g/cm , for five possible locations of the Moho-discontinuity, which 
is the boundary between the lower crust and the upper mantle. The bottom three rows 
give the real part of the first Fourier coefficient of the matter density distribution along 
the L = 1000, 1050, and 1100 km baselines. 

sponds the depth of 20 km at L ~ 500 km and 17 km at L ~ 800 km in Fig. 0, ±la (±3a) 
corresponds to the overall shifts by ±0.7 km (±2 km) of the Moho-discontinuity depth 
below the sea level. 

Because longer baselines go through the upper mantle for a longer distance, the 
average density grows with L; 2.854 g/cm 3 at L = 1000 km to 2.944, 2.980, 3.027, and 
3.051 g/cm at L(km) = 1050, 1100, 1150, and 1200, respectively. The average density at 
L = 1000 km is quite sensitive to the rise of the Moho-discontinuity because the baseline 
almost touches the mantle at L ~ 500 km for the mean estimate, as shown in Fig. [51 The 
average density grows from 2.854g/cm for the mean depth to 2.900g/cm , by 1.6%, if 
the Moho-discontinuity is only 0.7 km (la) higher than the present estimate. Although 
this is striking, the effect is significantly smaller than the 6% overall uncertainty in the 
conversion of the sound velocity to the matter density. We therefore neglect the effect 
in the present study, but further geophysical studies may be desired if L ~ 1000 km 
baseline is chosen for the far detector location. 

We also examine the impacts of the non-uniformity of the matter density distribution 
along the Tokai-to-Korea baseline. In Fig. [6J the black circles and white squares show, 
respectively, the real and imaginary parts of the Fourier coefficients of the matter density 
distribution along the baselines of L = 1000 km (a) and 1200km (b). As a reference, we 
show by the cross symbols the real part of the Fourier coefficients of the matter density 
distribution of PREM [27J. The PREM generally gives symmetric matter distribution 
about x = L/2, and hence does not give imaginary cofficients. When we compare them 
with the distribution for the Tokia-to-Kamioka baseline of Fig. HJ we find that Re(p~i) is 
now negative and its magnitude grows from 0.03 g/cm at L = 1000 km to 0.14 g/cm 
at L = 1200 km. We note here that Re(pi) at L = 1000 km depends strongly on the 
depth of the Moho-discontinuity, because whether the baseline goes through mantle or 
not affects the convexity of the matter density distribution significantly; see Fig. and 
discussions on the average matter density above. In order to show this sensitivity, we 
give Re(pi) values for L = 1000, 1050, and 1100 km in the bottom three rows of Tabled] 
It changes from — 0.03g/cm to — 0.08g/cm at L — 1000 km by more than a factor of 
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Figure 6: Fourier coefficients of the density distribution along the Tokai-to-Korea base- 
lines ; L = 1000 km (a), and 1200 km (b). Black circles and white squares show the real 
and imaginary parts of the Fourier coefficients, respectively, and the cross symbols show 
the real parts estimated from PREM [27J. 



two, if the Moho-discontinuity is only 0.7 km (la) shallower than our estimate of 20 km 
based on the measurement in Oki trough [2_8]. As for the L = 1050 km, the Re(pi) is 
almost insensitive to the depth of the Moho-discontinuity. Since the baseline of L — 1100 
km is more than 22 km below at Oki trough, the magnitude of Re(p~i) of L = 1100 km 
is insensitive to the exact depth of the Moho-discontinuity. 

Along the Tokai-to-Korea baseline, the convexity measure of — Re(pi)/po grows from 
~ 1% (~ 3% if the Moho-discontinuity is 3% shallower than our estimate) at L — 1000 
km to 5% at L — 1200 km. Their effects are hence non-negligible even with the 6% 
overall uncertainty in the average matter density. We also note that the negative sign of 
Im(pi) reflects the fact that the mantle is slightly closer to a far detector in Korea than 
J-PARC in Tokai. 



2.3 Comparison with PREM 

Before closing this section, we compare the matter density distribution of this work and 
those of the PREM [27], where the earth is a spherically symmetric ball such that the 
depth of the boundaries between adjacent layers are the same everywhere. In Fig. 
we show the cross section view of the T2KK experiment along the baselines according 
to PREM. Although the original PREM assumes that the sea covers whole of the earth 
down to 3 km from the sea level, we assume in this figure that none of the baselines go 
through the sea. 

There are several notable differences between the matter density distribution of this 
work shown in Fig. and that of PREM in Fig. [7J As for the T2K baseline, the Fossa 
Magna lies in the middle in Fig. |2l while it goes through only the upper crust in PREM. 
More importantly, the Conrad and the Moho- discontinuities are significantly shallower 
than those of PREM below Japan/East sea, whereas they are deeper than PREM below 
Japanese islands and Korean peninsula. Accordingly, although all but L = 1000 km 
Tokai-to-Korea baselines in Fig. Ogo through the upper mantle, the baselines of L £ 1100 
km are contained in the crust according to PREM. 
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Figure 7: The cross section view of the T2KK experiment along the baselines according 
to Preliminarily Reference Earth Model (PREM) [2j 
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PREM 


L (km) 


Po 
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Imp! 
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Po 


Repi 


Impx 


Rep 2 


295 


2.62 


+0.039 


-0.003 


-0.037 


2.60 











1000 


2.85 


-0.031 


-0.016 


-0.010 


2.71 


-0.098 





-0.040 


1050 


2.94 


-0.094 


-0.028 


+0.030 


2.76 


-0.094 





-0.016 


1100 


2.98 


-0.108 


-0.044 


+0.003 


2.78 


-0.090 





-0.030 


1150 


3.03 


-0.134 


-0.060 


-0.001 


2.93 


-0.201 





+0.040 


1200 


3.05 


-0.142 


-0.062 


-0.002 


3.00 


-0.220 





+0.009 



Table 2: Average and the Fourier coefficients of the matter density distribution along the 
Tokai-to-Kamioka baseline of Fig. [2] and the Tokai-to-Korea baselines in Fig. [5j which 
are compared with the predictions of PREM [27] . All the numbers are in units of g/cm . 

These differences are reflected in the average densities and in the Fourier coefficients 
of the matter distributions along the baselines. In Table [2J we summarize the average 
matter density (po) and a few Fourier coefficients along each baseline estimated from the 
matter density profile of this work, and those according to PREM, in units of g/cm . 

For the Tokai-to-Kamioka baseline, L = 295 km, despite the difference in the dis- 
tribution, the average matter density of 2.62g/cm turns out to be very similar to the 
upper crust density of 2.6 g/cm in PREM. This is because the matter density around 
Kamioka, being in the mountain area, is higher than the PREM upper crust density, 
which compensates for the low density in Fossa Magna. The presence of Fossa Magna 
around the middle of the baseline results in the positive value of Re(pi), reflecting the 
concave nature of the Tokai-to-Kamioka matter distribution. 

Along the Tokai-to-Korea baselines, all the average densities estimated from Fig. [5] 
are larger than those of PREM, because the crust density of 2.8g/cm in the region is 
higher than the world average of 2.6g/cm in PREM, and also because the traveling 
length in the upper mantle is longer than PREM, reflecting the shallower discontinuities 
below Japan/East sea in Fig. [5j Especially for L ~ 1050 km, the average matter density 
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of this work is larger than PREM by 6.1%. The difference reduces to 1.6% at L ~ 1200 
km. The negative values of Re(pi)/po represent the convexity of the matter density 
distribution. Their magnitudes are ~ 1% (~ 3% for 3% shallower Moho-discontinuity) 
at L = 1000 km and ~ 5% at L = 1200 km in this work, which are to be compared with 
~ 4% and ~ 7%, respectively, according to PREM. 

3 Neutrino oscillation in the matter 

In this section, we explain how we calculate the neutrino flavor oscillation probabilities 
with non-uniform matter distribution efficiently, and introduce analytic approximations 
that help our understandings on the matter effect. 

3.1 Exact evaluation of oscillation probabilities 

First, we show how we compute the oscillation probability exactly when the matter dis- 
tribution is approximated by step functions. In the earth, neutrinos interact coherently 
with electrons and nucleons via charged and neutral weak boson exchange, and these 
coherent interactions give rise to an additional potentials in the Hamiltonian. The po- 
tential from the neutral current interactions are flavor-blind, and it does not affect the 
neutrino flavor oscillation probabilities. Only the electron neutrinos, v e and z/ e , have 
coherent charged-current interactions with the electrons in the matter. This additional 
potential contributes only to the v e and v e phase, and hence affects the neutrino flavor 
oscillation probabilities [2"Tj . 

The location (x) dependent effective Hamiltonian of a neutrino propagating in the 
matter can hence be expressed as 



"<*> = h 




U \ m\ — m\ 



2 2 

mg _ m l 




(4) 



on the flavor space (i> e , z/ M , u T ) , after removing the components which are common to 
all the neutrino flavors. Here U is the lepton-flavor mixing matrix, the MNS matrix [TJ, 
defined as, 

Wa) = U ai \Vi) , (5) 

where u a represents the flavor eigenstate (u e , v^ and u T ) and ui represents the mass 
eigenstates with the mass m, (i = 1, 2, 3). The term a(x) in eq. PJ represents the matter 
effect; 



a\x) 



2V2G F En e (x) « 7.56 x lO^eV 2 ( -^3 J (t^T?) > (6) 

Vg/cm / VvjeV/ 



where Gp is the Fermi coupling constant, E denotes the neutrino energy, and n e (x) is 
the electron number density in the matter at a location x which is proportional to the 
matter density in the earth, p(x). 

Because of the location dependence of the matter density p, the Hamiltonian eq. (jlj) 
has explicit dependence on the location (x) of the neutrino along the baseline, or on the 
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elapsed (t = x/c) since the neutrino leaves J-PARC. This time dependence can be easily 
solved for the step-function approximation of the matter profile, since the Hamiltonian 
is time-independent within each layer of the region with a constant density[33]. The 
time evolution of the neutrino state can hence be obtained by the time-ordered product 
of the evaluations with time-independent Hamiltonian in each region; 



L n 

i/ Hdx =TT\ 



-iH k (x k -x k _ 1 ) 



(7) 



fe=i 



Here T denotes the time ordering, Hk is the time independent Hamiltonian in the region 
"Xk-i < x < Xk {k = 1 to n) along the baseline, Xk is the starting location of the region, 
where Xq = (km) is the location of the neutrino production point at J-PARC and 
x n = L (km) is the location of the target, SK or a far detector in Korea. 

Although it is straightforward to diagonalize the Hamiltonian eq. (J4J) in the region 
Xk-i < x < Xk where the matter effect term a(x) is a constant, we find that an ef- 
ficient numerical code can be developed by adopting a specific parametrization of the 
MNS matrix which allows us to decouple the dependence of the oscillation amplitude 
on the mixing angles 823 and the phase 5 from the matter effects [34]. The standard 
parametrization of the MNS matrix can be expressed as 



u = o 2 , P 13 P f o 



'12 



C23 S23 

— S23 C23 



p 



Cl3 



■«13 



Sl3 



Cl3 



Pt 



Cl2 «12 
-512 C12 



(8) 



by using a diagonal complex phase matrix P = diag (1, 1, e lS ), where Sij and c^- represent 
sin 6^- and cos 6ij, respectively. The Hamiltonian Hk of eq. (HI) with a(x) = a^ can be 
partially diagonalized as 



Hu 



— 23 P 
IE 



13 



13^12 



ms 



m\ 



1 2 

mj — m\ 



( 



OlOl + 



(ik 



P'ol 



2E 



:0 23 P 



O k 



(A2 — Aijfc 



(A3 — Ai)^ 



T k 



P ^23 + 2E , 



(9) 



where Xik (i = 1, 2, 3) are the eigenvalues of the matrix inside the big parenthesis and 
Ok is the ordinal diagonalizing matrix. Remarkable point of this parameterization is 
that the parameters #23 and 5 are manifestly independent of the matter effect [34] . The 
decomposition eq. (Q allows us to compose a fast and accurate numerical program for 
neutrino oscillation in the matter, since the non-trivial part of the Hamiltonian can be 
diagonalized by using only real numbers. 
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For the step-function like matter profile that gives eq. (J7|), the oscillation amplitude 
can be solved explicitly as 



S/3a(L) = {v p \ Te 4 Jo Hdx \v a ) 

n 

= (u p \ 23 PT J] [d k e- iAX ^*-^dl] pt0 2 T 3 \ Ua ) x (phase) . (10) 



fc=i 



Here 



AA fc = — 
IE 



(A2 — Aijfe 



(A 3 - Ai) fe j 



(11) 



gives the flavor non-universal part of the diagonal Hamiltonian in the region Xk~i < 
x < Xk where the matter density p and hence the matter effect term can be regarded 
as a constant a = a,k- Finally, we obtain the exact formula for the flavor transition 
probabilities, 

P u ^ u ,{L) = \S Pa {L)\ 2 . (12) 

It is explicitly seen in eq. (TTOj) that the dependence on #23 and 8 of the oscillation 
amplitude and the transition and survival probabilities are completely independent of 
the matter effect HE 



3.2 Perturbation formula 

The exact solution of the time evolution operator, eq. ( 1T2|) . is useful for numerical com- 
putation but it does not illuminate our understanding of the effect of the matter density 
distribution. In this subsection, we present a perturbation formula for the neutrino os- 
cillation probability with non-uniform matter profile. It helps us to understand not only 
the mechanism of determining the neutrino mass hierarchy pattern by using the average 
matter effect along the baselines, but also the effects of the matter density distribution 
along the baselines. 

Present experimental results tell that, among the three terms of (mass) 2 dimension 
in the Hamiltonian (J3J), \m\ — m\\ ~ 10~ 3 eV 2 is much larger than the other two terms, 
m 2 , — m\ ~ 10~ 4 eV 2 and a <■ 
splitting of the Hamiltonian 



10 eV for p ~ 3g/cm and E ~ 1 GeV. This suggests 
I]) into the main term H and the small term Hi. 



H(x) = H + Hi(x) 



(13a) 



H n 



IE 



:U 



HAx) = — 



1 
2E 



U 



1 1 
777.3 — m i 



777; 



777; 



C/t 



17+ 



a(x) 



(13b) 



13c) 
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When HqL is order unity, we can treat H\(x) as a small perturbation. We examine the 
second order approximation: 

Te -iJ o L Hdx ^ e -iH L _ • J d xe iHo(*-L) H f x \ e -iH x 

Jo 

-/ dx rdye iHo{x - L) H 1 (x)e iHo ^- x) H 1 (y)e- iHoy , (14) 

Jo Jo 

where H is location independent, and the location dependence of Hi(x) is due to the 
matter profile p(x), or a(x) along the baseline. 

The integral can be performed analytically by using the Fourier expansion of the 
matter profile along the baseline 



, N ^ ( 2nkx 

k=— oo \ / 



a + 2]T 
fc=i 



_ ,_ . {27ikx\ , . . {27ikx\ 

Re fofcj cos — - — + 1m (a fc ) sm 



L 



\ L ) 



(15) 



where 

* = 7. 56 x 10-,V (^) (^) (16) 

from eqs. d2J) and ([S])- We further divide /fi to the location independent part and the 
location dependent part 



H x (x) = H 1 + 6H 1 (x), 



( Oo 



oo 1 

8H x {x) = £ 



fc=i 



£ 



_, , . f2nkx\ _ , . . (2-Kkx\ 
Re (a fc ) cos I — — I + Im (a*) sm I — — I 



/ 1 



J V 



(17a) 
(17b) 



| . (17c) 




The v a — > up transition matrix elements can now be expressed as 

S(L) 0a = (up\ Te-^o^+x^* \ Va ) « (S (L) + S X (L) + S 2 (L))^ a , (18a) 

S (L)f ia = (uf i \e- iHoL \v a ) , (18b) 

S 1 (L) [Sa = -i (u p \ U L dxe iH ^^H 1 {x)e- iH A \v a ) , (18c) 
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S 2 {L) Pa = - (u fi \ (f L dx f X dye iH ^ x -^E x e iH ^ y -^E x e- iH ^\ \v a ) . (18d) 



Here we retain only the constant part H\ of the perturbation in the second-order term, 
eq. (118dl) . because the higher Fourier modes contribute negligibly in this order. The 
analytic expressions for the above amplitudes after integration are shown in the appendix. 
For the u^ survival mode, a — /3 — /i, we find 

S (L)^ = 1 + \U, 3 \ 2 (e' iAl3 - l) , (19a) 

2 a ;\tt \2\tt |2^0-^ ( -, , -iAi 



S X (L)^ = -i\U^\ 2 A 12 - ip^VW^^ (1 + e 



1 /_ ~ 2Re(a fc ) 



-2\U e3 \ 2 \U^\ 2 ^- 2 ao + E ; A 2 /A V* ( e_iAl3 " X ) ' ( 19b ) 
mi-mi\ £r{l-4k 2 (Tr/A 13 yj\ > 

S 2 (LU = -^f\U, 2 \ 2 + (a 2 ,aA 12 ) , (19c) 

where the phase Ay is 

m 2 — m 2 (m 2 — m 2 )\eV 2 } r , s 

In eq. fl!9cj) . we ignore the matter effect terms proportional to (a ) 2 because they are 
proportional to |£/ e 3| 2 , which is constrained to be smaller than 0.044 for \m 3 — m\\ = 
2.35 x 10~ 3 eV 2 by the CHOOZ experiment [Hj, and hence can be safely neglected in the 
second-order terms. The Fourier terms, a& in eq. ( 119bl) . come along with the factor of 
1/(1 — 4/c 2 (7r/A 13 ) 2 ), which suppresses contributions of the large k modes. 
The Up survival probability can now be estimated as 

P Vfi ^{L) » \S (L)^ + S^L)^ + S 2 (L)^\ 2 

= 1 - 4\U, 3 \ 2 (l - \U, 3 \ 2 ) (1 + A^sin 2 (^ + B^ + (a 2 , aA 12 ,A 2 2 ) , (21) 

where A^ and B^ are the correction terms to the amplitude and the oscillation phase, 
respectively; 



2 /_ >fi 2Re(a fc ) \ |E/ e3 | 2 (1 - 2|^ 3 | 2 ) A 2 2 W^\ 2 
ml- m\X l ° ^ 1 l-4F( 7 r/A 13 ) 2 y' 1 - \U^ 2 1 - \U^\ 



a L \U e3 \ 2 (l-2\U^\ 2 ) _ A 12 |Ey ^ 
AE 1-|C/^ 3 | 2 2 1-|C/ M 3| 



S =TF i 17TT2 5~i l777^- ( 22b ^ 
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In the limit of A^ = B^ = 0, the expression (12 ip reduces to the two-flavor oscillation 
probability in the vacuum. 

In order to show typical orders of magnitude of these correction terms, we evaluate 
them for the allowed range of the parameters [31 HI El El El E] : 



irrin — m 



2 1 

m 2 — m 1 



hO.ll 



-3„Ar2 



2.35^ IxHTW, 



7.50+°;^ x l(T 5 eV 2 



(23a) 
(23b) 



4|^ 3 | 2 (1 - l^sl 2 = sin 2 2^ ATM > 0.90 



(23c) 



4\U el U e2 \ 2 = sin 2 26 SOL = 0.852+°; 



026 ' 



(23d) 



4|[/ e3 | 2 (1 - \U e3 \ 2 ) = sm 2 26 RCT < 0.17 for |m 2 - m\\ 



2.35 x 10^eV 2 , (23e 



where in eqs. (I23cj) and (I23e[) . the 90% CL bounds are shown. From eqs. (123a[) . ( I23b|) 



and ( 23d ). we can set 



A12 = 0.032|A 



13 1 



(24) 



with 5 % accuracy, and the allowed range of 1 1/^3 1 2 and \U e3 \ 2 can be parametrized as 



\U^\ = sin 2 6 1 atm = 0.5 + O.lx. 
|^ e3 | 2 = sin 2 ^ RCT = 0.04y, 

with two real parameters x and y which can take values 



(25a) 
(25b) 



\x\ 



|^a| 2 -0.5 



0.1 



\tt r 
0.04 



< 1.6. 



(26a) 
(26b) 



in the 90% CL allowed regions of eqs. f)23cl) and (|23e|) . 

We now obtain the following expressions for N 1 and _B M around the oscillation max- 
imum, IA13I ~ 7r: 



N" ~ -0.0017 



7T L 



A) 



Ai 3 295km V3g/cni 



1 - 0.67- 



B^ ~ -0.033 [1 - 0.23^/y cos 5 



IA 



13 



MNS 



7T 



Re(p0y 


xy- 


/ A s 2 

- 0.0034 ( 13 J , 


(27a) 
(27b) 



From eq. f)27al) . we find that the main uncertainty in the oscillation amplitude comes 
from the sign of m 2 — m\ or A 13 in the term which is proportional to the matter density 
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and the baseline length. Its magnitude around the first oscillation maximum, IA13I ~ n, 
is constrained to be less than 0.26% for p = 2.62g/cm and Repi/po = 0.015 for the 
T2K experiment; see Table 2. This is about 4 times smaller than the proposed sensitivity 
of 1% for sin 2 26*atm in the T2K experiment [TH] . 

As for the corrections to the oscillation phase, the phase-shift B^ can be as large as 
0.04 for S = vr and y = 1.1 (sin 2 26> RCT = 0.17) at |Ai 3 | = vr. Although the T2K experi- 
ment is expected to measure the location of the oscillation maximum IA13/2 + B^\ = tt/2 
very precisely with the proposed error of 0.4% for |m 2 —to 2 ], the magnitude of \m 3 — m\\ 
should depend on the mass hierarchy which can differ from by as much as 3% [18J. 

For the v„ — > v e oscillation, Sq and Si for the v„ — > u e mode can be expressed as 



[S (L) + S 1 (L)} 



e/j. 



1 + 



rrin — m 



~ , . 2A 2 

a + 2^Re(<V 



fe=i 



Af 3 -4fc% 2 / 



[l-2\U e3 \ 2 ) 



-iJ2lm(a k ) 



4nkA 



13 



fc=l 



A 2 3 - 4fc 2 7T 2 



^3 (e~ iAl3 " 1 



+iU e3 U, 



e3^ M 3' 



2£ 



l-ltL-r U + e 



„-iAi 



iu e2 u; 2 A 12 . 



(28) 



The remarkable point is that the contribution of the Im(ofc) terms to the transition 
probabilities \S + 5*1 1 2 is highly suppressed, because they cannot interfere with the 
leading So term and hence they contribute to the probability only as terms of order 
a or A 12 . The oscillation probabilities are hence insensitive to the asymmetry in the 
matter density distribution when the sub-leading phase, A 12 is small. We confirm this 
observation quantitatively in the numerical analysis and neglect all terms of Im (a^) in 
the following. 

The transition probability can then be expressed as 



Pu^uXL) « |So(L) e/ , + S 1 (L) e/ ,| 2 + 2Re(S (L) e ^ 2 *(L) 



(-11 j 



4|^3r|f/e3| 2 {(l + ^ e )sin 2 ^) +J B e sin(A 13 )} + C e , 



(29) 



where A e and B e , and C e are the correction terms; 
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(30b) 
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A 12 Re U e2 U; 2 U* e3 U^ ( - ao \ tflva u e2 u; 2 u*jj, 



1 + 



'fi3 



1 + ao 



I^H^I 2 V ml -m\ 4 |^ 3 | 2 |^e3| 2 V* m\-m\ 



■Tlrll ,/ ,/ ,, , f^l 1 ' 
2£ 



C e = A 2 12 \U e2 \ 2 \U ll2 \ 2 -2A 12 ^Re hT rr* rr* rr " , /^Yir; ,2,,, ,2 



Ue2U; 2 U* e3 U ll3 \ + ^j \U e3 \ 2 \U,s\ 2 ■ (30c) 



When we substitute the present experimental constraints of eqs. ( 123]) - ( 125]) . they can be 
expressed as 

* L Po h nno ^f■, , ^ 2Re(p!/po) \ 



/4 e ~ 0.11- ^ (1 - 0.08y) [l + y 
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vr \ ' A 13 295km 3g/cni 

sin<5 M NS /Ai 3 \ 2 / 7T L p 



1 : 1-7 , A , onc1 ,„.„ 7 3 J , (31b) 

A 13 \" nn#VIO _ . IA13I L , „„ { L \ 



y/y V 7r / V |A 13 | 295km 3g/crn 

2 I A I r / r \ 2 



C^aOOU^j - .0013^oos* MNS ^— + 0.00036^— j ,(31c) 

where we set x = df/^] 2 = 0.5) in the correction factors above for the sake of brevity. 
Finite values of x in the allowed values of |x| < 1.6 do not affect the following discussions 
significantly [TH]. When the magnitudes of A e and B e are both much smaller than unity, 
these two terms can be factorized as eq. ( 1211) for the z/ M — > v^ probability. Therefore 
A e affects the amplitude of the oscillation and B e gives the oscillation phase shift from 
A13/2, except for the small term C e , which does not vanish in the limit of \U e ^\ 2 — > 0. 
It is clear from the above parametrization that the magnitude of A e and B e can be as 
large as 0.5 or bigger for L £ 1000 km, and we find the expression eq. (1291) more accurate 
than the factorized form. 

The first term in A e in eq. (I31al) gives the matter effect, whose sign depends on the 
mass hierarchy pattern. When the hierarchy is normal (inverted), the magnitude of the 
v \l -^ v e transition probability is enhanced (suppressed) by 10% at Kamioka, and by 30 
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% at L ~ 1000 km in Korea around the oscillation maximum, IA13I ~ n. If we measure 
the maximum of the oscillation probability ~ 4|[/ /i 3| 2 |[/ e 3| 2 (l + Ap Cak ) at two locations, 
we should expect 



A^ k (L ~ 1000km) - A^(L = 295km) ~ \ ™ £ ™ Z^'™ ^ (32) 



-0.2 for the normal hierarchy , 
-0.2 for the inverted hierarchy , 



within the allowed range of the model parameters. Therefore, simply by comparing the 
magnitude of the oscillation maximum between SK and a detector at L ~ 1000 km, we 
can determine the neutrino mass hierarchy [T5]. Once the sign of A 13 is fixed, sin<5 can 
be determined via the third term of eq. (13 lap . 

Fourier modes of the matter distribution affect the magnitude of the first term in 
A e , and hence affects the sin 6 measurement. Since the factor 1/(1 — 4k 2 (iT /A13) 2 ) is 
negative around the first oscillation maximum, negative Ke(pk) enhances the average 
matter effect, and positive Re(pfc) suppresses it. Figs. EJEJ and Table [7] show that Re(pi) 
is positive for the Tokai-to-Kamioka baseline, while it is negative for the Tokai-to-Korea 
baselines; see Table |2j Hence the first Fourier mode reduces the average matter effect 
very slightly by ~ 1% for the T2K experiment, while it enhances the effect by about 
1 ~ 5% for the Tokai-to-Korea baselines at around |A 13 | ~ n. Contributions from the 
higher Fourier modes are highly suppressed. 

The term 5 e sinAi3 in eq. ( 129]) is also sensitive to the sign of A13, and hence the 
mass hierarchy. Although its sign depends on the magnitude of cos5/y/y = cos5/\U e3 \ 
in the second term of eq. (131b|) . the difference in its magnitude at two locations satisfy 



B e (L ~ 1000km) - B e (L ~ 295km) 0.2 , (33) 

within the allowed range of the model parameters. If we define the phase of the oscillation 
maximum by 

|A^ ak (L) + 25 e (L)| =tt, (34) 

we should expect 



^peak^ ^ 1000km) - A pcak (LKL = 295km) ~ { + °' 4 for the normal hierarchy, , . 
13 13 ) —0.4 for the inverted hierarchy . 



This implies that the mass hierarchy can be determined also by measuring the location 
of the v^ — y v e oscillation maximum with the accuracy better than around 10%. 

It should also be noted that the terms proportional to sin S and cos S in eqs. (13 lap and 
( I31bp are proportional to l/^/y oc l/|£/ e 3|. Because the v^ — > v e oscillation probability 
is proportional to y oc \U e s\ 2 in eq. ( 1291) . the statistical error of its measurement is 
proportional to ^fy. The 1/y 7 ^ factor in front of the sin S and cos 5 terms cancel this ^Jy 
factor, and hence the error of the S measurements does not depend on the magnitude of 
|f/ e 3| 2 , as long as the neutrino mass hierarchy is determined [T7J [18] . 

From the above observations, we understand quantitatively the reason why the mass 
hierarchy and both sin 6 and cos 8 can be measured by observing v^ — y v e oscillation 
around the first oscillation maximum at two vastly different baseline lengths. The T2KK 
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proposal of Refs. [T7J [18] showed that this goal can be achieved by using the existing 
detector SK at L = 295 km and only one new detector at a specific location in the 
Korean east coast where the beam can be observed at an off-axis angle below 1°, so that 
the neutrino flux at the first oscillation maximum (~ a few GeV) is significant. 

In this report, we investigate further the impacts of using v^ beam in addition to v^ 
beam. The oscillation probabilities Piv^ — > v^) and Piy^ — > v e ) are obtained from the 
expressions for P(z/ M — y u^) and Piy^ -y z/ e ), respectively, by reversing the sign of the 
matter effect term (p — y —p) and that of 5. Because these sign changes occur indepen- 
dently of the L-dependence of the oscillation probabilities, we may expect improvements 
in the physics potential of the T2KK experiment for a given number of POT. 

3.3 Hybrid method 

In this subsection, we introduce a very simple approximation formula for the oscillation 
probabilities that can be used for fast simulation purposes. The analytic formula of 
the previous subsection are useful in understanding qualitatively the dependences of the 
oscillation probabilities on the model parameters and the baseline lengths, but they are 
not suited for quantitative analysis presented below. It is because the approximation of 
keeping only the first and second orders of the matter effect terms pi and the non- leading 
oscillation phase A12 is not excellent at far distances (L ~ 1000 km), as can be inferred 
from the magnitude of the correction terms A e and B e in eq. ( )30l) . 

On the other hand, the average and Fourier coefficients of the matter density profile 
shown in Table [2] for Tokai-to-Korea baselines tell that Re(pi)/po is rather small, 1.5% 
for T2K and 1 ~ 5% for Tokai-to-Korea baselines (L = 1000 ~ 1200km). Higher 
Fourier modes are expected to contribute even less. Therefore, we expect that the 
perturbation in terms of Re(ai) to the oscillation probabilities for the uniform matter 
density, p(x) = po, can give a good approximation. The v^ — y v^ survival and v^ —¥ v e 
transition probabilities in such approximation are expressed as 

PF% (So, Re(a l} ) = P~% (ao) + ^j ^^ ^ - - ^ sin 2 (*£) ,(36a) 
v^v^\ u, v v)) „„-►,,„ v uj m 2_ m 2 i-4(vr/A 13 ) 2 V 2 J ' v ; 



phybrid /- p /- X, pcxact (- S ^ Iggggl) |^e3^3 ? (1 ~ 2|^ e3 | 2 ) . 2 /A 



(36b) 



respectively. We call the approximation hybrid, being a combination of the exact oscil- 
lation probabilities for a constant matter density and the first order analytic correction 
term which is proportional to the real part of the first Fourier coefficient of the matter 
density distribution along the baseline. 

In eq. ( I36a[) . the second term is suppressed because |1 — 2|[/ At3 | 2 | = 0.2x < 0.32, 
see eqs. (125]) and (f26j) . We find that eq. (I36aj) gives an excellent approximation to the 
v p, ~^ v n survival probability in the whole range of the parameter space explored in this 
report. Equation (136bjl also gives an excellent approximation to the v^ — y v e oscillation 
probability around the oscillation maxima. However it does not always reproduce accu- 
rately the probability around the oscillation minima at which the first term gets very 
small and our first order correction term can be off significantly due to the phase-shift. 
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We find that the difference between the exact formula and the hybrid method around the 
oscillation minimum can be as large as 10 ~ 30%. Because the number of event is also 
small around the minimum region, this difference is not significant for the A% 2 analysis 
in most of the parameter space. We find that the minimum A% 2 values are obtained 
rather accurately when the exact formula are replaced by the approximation of eq. ( 136|) . 

4 Analysis method 

Before we present the results of our numerical calculation, we would like to explain our 
treatment of the signals and backgrounds, and the analysis method. 

In order to compare our results with those of the previous works, we adopt the same 
setting of Ref. [18J. There, the detectors at both Kamioka and Korea are assumed 
to have excellent detection and kinematical reconstruction capabilities for v^ and v e 
Charged Current Quasi-Elastic (CCQE) events within the fiducial volumes of the 22.5 
kt at Kamioka (SK) and 100 kt at Korea. The energy bin width, 5E, in the analysis 
is chosen as 200 MeV for E u > 400 MeV, which accounts partially for kinematical 
reconstruction errors 2 . 

The v^ — > v a oscillation signals in the z-th energy bin, E\, = 200MeV x (i + 1) < 
E v < E l u + 5E U , at each site are calculated as 

Nl D (i/ M ) = M D N A / . $^(E, L) P v ^ Va {E) v C a CQE {E) dE , (37) 

J El 

where the suffix D = SK or Kr denotes the detector, and P Ufi ^ Ua is the neutrino oscillation 
probability, which is calculated by the exact formula, eq. (JTj). Msk = 22.5 kt and 
Mxr = 100 kt are the detector mass within the fiducial volume, Na = 6.017 x 10 23 is the 
Avogadro number, $„ is the v^ flux from J-PARC [35] which is proportional to 1/L 2 , 
and a^ CQE and erf CQE are, respectively, u^ and u e CCQE cross sections per nucleon in 
water [SB] . 

As for the background, we only consider the contribution from the secondary neutrino 
fluxes of the u^ primary beam, such as u e , z/ M , u e which are calculated as in eq. (1371) where 
$ Vli (E) and P v ^ Va are replaced by $ Vf) (E) and P Uli ^u a , respectively, for up = v e , v e , u^, 
and u a = u e or u e for a = e, and u^ or u^ for a = //. The data of all the primary and the 
secondary fluxes as well as the CCQE cross sections used in our analysis can be obtained 
from the website [22] • After summing up these background events, the numbers of e-like 
and /i-like events in the i-th bin are calculated as 

<d = K,d(^) + N «% G > (a = e, /*; D = SK, Kr) . (38) 

When the primary u„ beam is replaced by the primary u„ beam, the signals are obtained 



U e ,V P ,V, 



/-'• 



from Pp — > u a (a — e, (x) while the backgrounds are from up —¥ u a with up 

As a quantitative indicator of the capability of the proposed experiments to measure 
the neutrino model parameters, we introduce a Ax 2 function as follows: 

A X 2 = Xstat + Xsys + X?ara ■ (39) 



2 See e.g. Ref. 20 for a more realistic treatment, where smearing effects due to nucleonic Fermi 
motion, resonance production and finite detector resolutions have been studied. 
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The first term, Xstat> accounts for the parameter dependence of the fit to the CCQE 
events, 



xL t = £ E { I - e ' W- e ' 1 + 1 v ^== M ' 1 }i ( 4 o) 

D=SK,Kr i 



where the summation is over all bins from 0.4 GeV to 5.0 GeV for N^ at both sites, 
0.4 GeV to 1.2 GeV for V e)SK , and 0.4 GeV to 2.8 GeV for N e , Kl . Here N l eD and A^ D 
are the calculated number of e-like and /z-like events in the i-th bin in each detector 
(D = SK, Kr), and their square-roots give the statistical error. The neutrino oscillation 
probabilities P u ^ Ua (E) in eq. (J3TJ) are calculated for the fallowing input parameters: 

\m\ -m\\ = 2.35 x 10" 3 eV 2 , (41a) 

m 2 2 - m\ = 7.5 x 10" 5 eV 2 , (41b) 

|^ 3 | 2 = sin 2 #ATM = 0.5, (41c) 

MU e iU e2 \ 2 = sin 2 2#sol = 0.85 , (41d) 

for both hierarchies (mj — m\ > or m 2 — m\ < 0), and for various values of sin 2 2#rct 
and 5. For the matter density profile along the baselines, we use the mean matter density 
of each regions as given in Figs. [2] and [5] 

/ . \ fit / . \ fit 

The numbers, (^d) and IN' l lD ) are then calculated by allowing the model pa- 
rameters to vary freely and by allowing for systematic errors. We consider the following 
systematic errors in this analysis. We assign 6% uncertainty to the overall matter den- 
sity along each baseline, Tokai-to-Kamioka and Tokai-to-Korea, because the ambiguity 
in the density-velocity conversion scale dominates the uncertainty of the average matter 
density: 

PW™ = / P SK P SK (x) , p(*)™ = /f p Kr (x) , /f, ff = 1 ± 0.06 . (42) 

Although we expect positive correlation between the scale factors of fp and f^ v , we 
treat them independently as a conservative estimate. We assign 3% uncertainty in the 
normalizations of each beam flux: 

$%{E, L) = f Ua $ a (E, L), f Va = 1 ± 0.03 (for u a = Vfl , v^ z/ e , z/ e ) . (43) 

Here also, we ignore possible correlations among the errors of primary and secondary 
beam fluxes. For the CCQE cross sections of neutrinos and anti-neutrinos, we assign 
common 3% error for v^ and v e events; 

< CQE ' fit (£) = /< CQE (£), ^ c ^ fit (^) = /^ C e CQE (^), /, = 1±0.03, (44) 
and another common 3% error for v^ and v e events; 

< CQE ' fit (^) = A< CQE (£), ^^\E) = /^ C e CQE (^) , A=l±0.03. (45) 
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This is equivalent of assuming the e-fi universality in the cross section uncertainly, while 
neglecting correlation between v\ and v\ cross section errors. For the fiducial volume of 
SK and a far detector in Korea, we assign 3% error each: 



A*d = /dM d , f D = 1 ± 0.03 (for D = SK, Kr) . 



(46) 



Summing up, we introduce 10 normalization errors, when calculating (N l ) , and x 2 ys is 
expressed as 




+ E 

D=SK, Kr 



(47) 



Finally, XjLra accounts for external constraints on the model parameters: 



X 



para 



■ 2 - •'■■ lli ' - 7.5 x 10- 5 eV 2 \ 2 /sin 2 20&T - 0.85 N 



\m% — m 



0.2 x 10" 5 



+ 



0.025 



+ 



/sin 2 2^* -sin 2 20 W 






'K('.r " " ■">" -"[;CT_ ] 



0.01 






(48) 



The first two terms are essentially the present constraints from the KamLAND experi- 
ment [6], summarized in eq. ( 123]) . In the last term, we assume that the planned future 
reactor experiments [TQl El EG] should measure sin 2 26'rct with an uncertainty of 0.01, 
during the T2KK experimental period. We do not impose the present constraints on 
l m i — m il an d sin 2 26* atm given in eq. (1231) . since the sensitivity of the T2KK experi- 
ments supersedes them. 

In total, Ax 2 depends on six model parameters and ten normalization factors. The 
minimum of Ax 2 is then found in this 16 dimensional parameter space. 



5 The earth matter effects and the mass hierarchy 

In this section, we discuss the relation between the earth matter effects and the capability 
of determining the neutrino mass hierarchy pattern in the T2KK experiment. We also 
examine the dependence of the capability on the location of a far detector in Korea. 



u^ and 



5.1 Determining the mass hierarchy 

The analytic expressions presented in the previous section show that v^ — >■ 
Vy. ~^ ^ijl oscillations are insensitive to the neutrino mass hierarchy, while v^ — > v e 
and v^ —¥ v e oscillation probabilities depend significantly on the mass hierarchy pattern 
through the terms which are proportional to the baseline length L and the average matter 



26 



5 


0° 


90° 


180° 


270° 


previous results [18] 


14.6 


23.0 


54.7 


42.9 


changes in input model parameters 


14.7 


25.1 


53.1 


48.7 


overall error of the density (3% — > 6%) 


14.4 


24.7 


51.2 


47.0 


average density SK(2.8g/cm — >■ 2.6g/cm ) 


14.6 


24.9 


51.8 


47.0 


average density Kr (3.00g/cm — » 2.98g/cm ) 


14.2 


24.4 


51.0 


47.0 


non-uniform matter density profile 


15.7 


26.8 


54.2 


49.1 



Table 3: The changes in minimum A% 2 values at each modification on the matter effect 
treatment from the previous work [18J. Signals are generated for sin 2 2#rct — 0.10 and 
the four 5 values by assuming the normal hierarchy and fit is performed by assuming the 
inverted hierarchy. The combination of off- axis angles are 3.0° at SK and 1.0° at Korea, 
and the Tokai-to-Korea baseline length is 1100 km. The results are for the neutrino 
beam with 5 x 10 21 POT. 



density po- Two independent observables are sensitive to the mass hierarchy. According 
to the estimates of eqs. (132]) and (|35l) . the magnitude of the first oscillation maximum at 
L ~ 1000 km should be about 20% bigger (smaller) than that of L = 295 km when the 
mass hierarchy is normal (inverted). Likewise the location of the phase where the first 
peak occurs at L ~ 1000 km is about 25° bigger (smaller) than that at L — 295 km. All 
the above predictions should be reversed for v^ — > v e oscillation. It is remarkable that 
none of the other model parameters contribute significantly to the differences eqs. (1321) 
and (1351) . since they tend to cancel between the two measurements of the same oscillation 
probabilities. 

Since the matter effect is proportional to the average matter density po along each 
baseline with corrections from higher Fourier modes, we expect that capability of con- 
straining the mass hierarchy pattern in the T2KK experiment depends on the earth 
matter profile along the baseline. 

In order to examine the impacts of the non-uniform matter profile along the baselines, 
let us compare our result with those obtained in Ref. [IB], where the matter profiles 
along the baselines are assumed to be uniform with densities 2.8g/cm 3 for T2K and 
3.0g/cm 3 for Tokai-to-Korea. The errors of these densities have been assumed to be 
3% [T71 HH US]- We summarize changes in A% 2 values for each modification step by 
step in Table. EJ Here, the minimum A% 2 is obtained from the data generated by 
assuming the normal hierarchy (m 2 — m 2 > 0) while the fit is performed by assuming 
the inverted hierarchy (m 2 — m 2 < 0). The input parameters are sin 2 2#rct — 0.10 and 
5 = n x 90°(n = 0, 1, 2, 3). The combination of off-axis beam is 3.0° at SK and 1.0° 
in Korea at L = 1100 km. 

The first row of Table [3] gives input 5 values, and the second row gives the results 
of Ref. [TS] which are reproduced in our analysis by using the same input parameters. 
In the third row, exactly the same computation is repeated when the input parameter 
values are updated according to more recent experimental data, summarized in eq. (123]) . 
It is remarkable that the minimum A% 2 values do not change much even though the 
errors of m 2 , — m 2 and sin 2 2#sol have been reduced by a factor of three since 2006 when 
the analysis of Ref. [18] were performed. 

The first modification in the matter effect is the magnitude of the uncertainty in the 
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matter density, which changes from 3% in Ref. jTSj to 6% in this report according to 
discussions in section 2; see Fig. |3j Although the matter effect is the key to determine 
the mass hierarchy pattern in the T2KK experiment, the larger error of the density 
does not affect A% 2 significantly. This is because the statistical error is dominant in 
the present setting. We note in passing that we assign 6% error independently for T2K 
(/p K ) and Tokai-to-Korea (/? r ) baselines, even though the dominant systematic error 
in the conversion of sound velocity to the matter density is common. This makes our 
estimates conservative, since the typical size of the independent error is ~ 3% [16J. 

The second and the third modifications are the values of the average matter densities. 
The reduction of the average matter density for the T2K baseline due to Fossa Magna 
improves the mass hierarchy discriminating power since the differences eqs. ( 132]) and 
(|35|) increase. A slightly smaller value of the average matter density at L ~ 1100 km 
decreases A% 2 : The latter effect is non- negligible because the differences are proportional 
to 1100 x p (L = 1100km) - 295 x p (T2K). 

Final modification is the non-uniformity of the matter density distribution. It is 
remarkable that the minimum Ax 2 increases significantly for all four 5 cases, as conse- 
quences of non-uniformity in the matter density distributions. The effect is particularly 
significant for T2KK, because Re(pi) is positive for the T2K baseline due to Fossa Magna; 
see Fig. M and Table d 

In Fig. [HJ we show the minimum values of A% 2 for the T2KK two detector experiment 
as functions of the matter density distributions which are approximated by the distri- 
bution with the average density (average), including the real-part of the first Fourier 
mode (Repi), further including the imaginary part (Impi), and also with the second 
modes (p 2 ). The results are then compared with Ax^in calculated by using the exact 
(step-function-like) distribution. As a reference, A^ 2 ^ calculated by the hybrid ap- 
proximation of the oscillation probability, eq. ( )36l) . is also shown (hybrid). The input 
hierarchy is normal in the left figure while it is inverted in the right figure, where in both 
cases the opposite hierarchy is assumed in the fit. The input sin 2 2#rct is 0-10 f° r solid 
red blobs, 0.06 for open blue symbols, whereas 0.02 for solid green blobs, whereas the 
input CP phase are 0° for squares, 90° for upper triangles, 180° for circle, and 270° for 
upside-down triangles. We show the results for a far detector at L = 1100 km baseline 
observing 1.0° OAB, which has — Re(pi)/po = 0.036. Fig. [8] shows that the minimum 
A% 2 values are reproduced accurately for all the cases just by accounting for the Re(pi) 
component. We conclude that higher Fourier modes including Im(pi) do not affect the 
minimum A% 2 significantly in the T2KK experiment. 

The approximation of keeping only the linear term in Re (pi) in the oscillation prob- 
abilities, the hybrid method of section 3.3, works rather well at most cases. However, 
in Fig. Mb), the hybrid method gives slightly high minimum A% 2 for 5 = 180° and 90° 
at sin 2 2# RCT = 0.1, and for 5 = 180° at sin 2 26 I RCT = 0.06. In these cases, the hybrid 
method overestimates A% 2 in bins around the oscillation minimum where the approxi- 
mation is poor; see discussion in section 3.3. In more realistic studies including neutral 
current 7r° background, see e.g. [20], those bins around the oscillation minima play less 
significant role due to the background dominance, and we expect that the hybrid method 
gives more reliable results. 
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Figure 8: Minimum A% 2 of the T2KK two detector experiment when the matter distri- 
bution along the baselines are approximated by their average, and up to first and second 
Fourier modes. They are compared with the results of the exact (step-function-like) mat- 
ter distribution and those of the hybrid approximation to the oscillation probabilities. 
The input hierarchy is normal (inverted) in the plot a (b) where the opposite hierarchy 
is assumed in the fit. The input sin 2 26'rct are 0.10 for solid red, 0.06 for open blue, and 
0.02 for solid green symbols, and the input CP phase are 0° for squares, 90° for upper 
triangles, 180° for circle, and 270° for upside-down triangles. The results are for 3.0° 
OAB at SK, and for a 100 kt far detector at L = 1100 km observing 1.0° OAB, with the 
neutrino flux from 5 x 10 21 POT at J-PARC. 



5.2 Optimal T2KK setting for determining the mass hierarchy 

Since the matter profile along the Tokai-to-Korea baseline depends on the baseline length 
as shown in Fig. [5j we re-examine the location dependence of the far detector capability 
to determine the neutrino mass hierarchy pattern. We show the result of our numerical 
calculation in Fig. [9J which shows the minimum A% 2 to reject the wrong mass hierarchy 
for various combination of the off-axis angle and the baseline length of Tokai-to-Korea 
baseline, when the 3.0° OAB reaches SK. The left figure (a) is for the normal hierarchy 
and the right figure (b) is for the inverted hierarchy. The input parameters are \m\ — 



m 



2.35 x 10~ W 



rrin 



iin 



7.5 x 10-W^, sin^ATM = 0.5, sin^soL = 0.85, 
sin 2 26'rct = 0.10 and 5mns = 0°, and A% 2 is calculated for a far detector of 100 kt 
fiducial volume and 5 x 10 21 POT at J-PARC. It should be noted here that the 0.5° 
OAB reaches Korean peninsula only for L ;$ 1030 km; see the contour plot Fig. 1 of 
Ref. [17J. We therefore have only the L = 1000 km point for 0.5°, three points for up to 
L = 1100 km for 1.0°, all five points for higher off-axis angles in Fig. [9j 

It is clear from the figure that the combination of OAB 0.5° at L = 1000 km is most 
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(b) hierarchy: inverted 
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Figure 9: The minimum A% 2 to exclude the wrong neutrino mass hierarchy as functions 
of the off- axis angle and the length of the far detector in Korea, when the 3.0° OAB 
reaches SK. The left figure (a) for the normal hierarchy, and the right figure (b) is for 
the inverted hierarchy. The input parameters are the |m 2 — m 2 \ = 2.35 x 10~ 3 eV 2 
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8.3 x 10~ 3 eV 2 , sin 2 # 
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0.85, 



sin^2# RCT = 0.10 and 



^mns = 0°. Ax 2 is calculated for a far detector for with 100 kt fiducial volume and for 
the neutrino flux with 5 x 10 21 POT at J-PARC. 



powerful to determine the mass hierarchy pattern, for both the normal and inverted 
hierarchy cases. This is because the 0.5° OAB has the strongest flux around the first 
oscillation maximum [18]. We also find that the 1.0° OAB has significant sensitivity to 
the mass hierarchy especially at longer baseline lengths. The value of the A% 2 for 1.0° 
OAB at L — 1100 km is 15.7 (15.8) for the normal (inverted) hierarchy, respectively, as 
compared to 21.0 (18.2) for 0.5° OAB at L — 1000 km. On the other hand, the sensitivity 
of T2KK experiment on the mass hierarchy decreases significantly for the off-axis angle 
greater than about 1.5°. 

We also examin the case when the 2.5° OAB reaches SK, whose results are shown in 
Fig. [TOj It is worth noting here that no location inside Korean peninsula can observe the 
neutrino beam at an off-axis angle ~ 0.5° : See the off-axis angle contour plot of Fig 1 
in Ref. [TTJ. The physical location of a far detector that detects 1.0° OAB in Fig. [TUl is 
the same as that of observing 0.5° OAB in Fig. [9j In general, when we decreases the 
off-axis angle at SK by 0.5°, that of a far detector in Korea increases by ~ 0.5° as can be 
seen from the schematic picture of Fig. [TJ Therefore, by comparing the minimum A% 2 
values of the corresponding points in Fig. M and (TUl, we can estimate the off-axis angle 
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Figure 10: The same as Fig. 9 but for the case when 2.5° OAB reaches SK. 

dependence of the T2KK experiment for each location of the far detector in Korea. We 
generally find that the minimum A% 2 value decreases by about 40% when the off-axis 
angle at SK is decreased from 3.0° to 2.5°. It is therefore important that the off-axis 
angle at SK should be increased, i.e. the beam center should be oriented deeper into 
underground at J- PARC. 

6 Earth matter effect for the anti-neutrino beam 

In this section, we study the impacts of using the anti-neutrino beam in addition to the 
neutrino beam in the T2KK neutrino oscillation experiment. 

Potential usefulness of having both neutrino and anti-neutrino beam run in the two 
detector experiments like T2KK can be understood as follows. The leading contributions 
to the Vp — > u e and v^ — > v e oscillation shifts in the amplitude and the phase are 
expressed, respectively, as 



A e ~ ± 0.11 



7T 



L 



Po 



A 13 295km 3g/cm 



, 0.29*^ 



Vy J 



B e ~ ^0.08 



L po /cos S 

' +0.14(^^-0.11 



295km 3g/cm 



Vv 



(49a) 
(49b) 
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where the upper signs are for v^ — » v e and the lower signs are for v^ — > v e oscillations, and 
only the leading terms in eq. ( l3Tj) are kept. Because the matter effect terms proportional 
to LpQ contribute with opposite sign for v^ and v^ oscillations, we can generally expect 
higher sensitivity to the mass hierarchy by comparing the oscillation probabilities at large 
L. In addition, since only the matter effect term changes the sign in the phase-shift term 
B e of eq. (149bj) . possible correlation between cos S and the mass hierarchy pattern can 
be resolved, and we can expect improvements in the hierarchy resolving power. On the 
other hand, since the v^ flux is generally lower than the v^ flux for proton synchrotron 
based super beams, and also because the background level from the secondary v e flux is 
higher than the corresponding v e background for the v^ —¥ v e signal, we should evaluate 
pros and cons quantitatively. 

In order to estimate the merit of using both v^ and v^ beams at T2KK, we repeat 
the Ax 2 analysis by assuming 2.5 x 10 21 POT (2.5 years for the nominal T2K beam 
intensity [13]) each for v^ and v^ beams, so that the total number of POT remains the 
same. We allow independent flux normalization errors of 3% for the primary v^ beam, 
and the secondary z/ e , v^ and v e beams. The systematic error term of x 2 ys i n ec L- (|39|) 
now depends on 14 normalization factors; 



X 2 ^ = Xs 2 j 47) + E 1^1 . (50) 



where Xsys 47 ^ * s defined in eq. (147j) . and the four f a terms measure the flux normalization 
uncertainty of the primary (p^) and the secondary (u^, z/ e , u e ) fluxes in the anti-neutrino 
beam from J-PARC [H]. 

The results are shown in Fig. El for 3.0° OAB at SK, and Fig. [TJfor the 2.5° OAB at 
SK. It is remarkable that the splitting of the total neutrino flux into half v^ and half v^ 
results in significant improvements in the mass hierarchy resolving power of the T2KK 
experiment. For the combination of 3.0° OAB at SK and 0.5° OAB at L — 1000 km, the 
increase in Ax 2 ^ reads 21.0 to 33.4 for the normal hierarchy while 18.2 to 30.7 for the 
inverted hierarchy. 

It is even more striking that the hierarchy resolving power at the 3 locations of the 
far detector at 1.0° OAB is as high as that of L = 1000 km at 0.5 OAB in Fig. [TQ 
At L = 1100 km (blue open triangles), Ax^ increases from 15.7 and 15.8 in Fig. EH 
respectively, to 33.9 and 40.7 in Fig. [11] for the normal (a) and inverted (b) hierarchy 
case. We find that this is partly because of our choice of the CP phase, 5 = 0°. For 
cos S = 1 and sin 2 26 l RCT = 0.1, the phase-shift term B e in eq. f!31bj) becomes as large 
as 0.4 for the v^ — > v e oscillation, significantly shifting the peak location of eq. (34) 
to larger |A 13 |, or smaller E„, for the inverted hierarchy (A 13 < 0). This shift of the 
oscillation maximum to lower E v compensates for the slightly softer spectrum of the 1.0° 
flux as compared to the 0.5° OAB flux. We confirm the expected S dependence of this 
observation by repeating the analysis for S = 180°, where the increase in A^ 2 ^ is not 
as dramatic as in the 6 = 0° case. General survey of the physics capability of the T2KK 
experiment over the whole parameter space of the three neutrino model is beyond the 
scope of this report, and will be reported elsewhere. 

We show in Fig. [12] the corresponding Ax^^ values when the beam center is oriented 
upward by 0.5° so that the 2.5° OAB reaches SK, and hence the far detector in Korea 
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Figure 11: The same as Fig. |9] but when both neutrino and anti-neutrino beams are 
used. Ax 2 is calculated with 2.5 x 10 21 POT for each beam so that their sum agrees 
with the neutrino flux assumed in Figs. |9] and [T0J 

observes the beam at 0.5° larger off- axis angles. When we compare the A^ 2 ^ values 
in Figs. [IT] and d2] at the same far detector locations, we generally find smaller values 
in Fig. H2J which are consequences of the softer neutrino beam spectra at larger off-axis 
angles. As in the neutrino beam only case, it is desirable to make the off-axis angle as 
large as possible at SK, so that the far detector in Korea can observe the same beam at 
a smaller off- axis angle. 

7 Summary and discussion 

In this paper, we study the earth matter effects in the T2KK experiment by using recent 
geophysical measurements [28j EH [30j |3TJ [32] . 

The mean matter density po along the Tokai-to-Kamioka baseline is found to be 
2.6g/cm . The presence of Fossa Magna along the T2K baseline makes the matter density 
distribution concave, or Re(pi) > for the first coefficient of the Fourier expansion 
of the density distribution along the baseline, which is opposite from what is naively 
expected from the spherically symmetric model of the earth matter distribution such as 
PREM [27|. We find that both the reduction of the average matter density po and the 
positive Re(pi) contribute positively to the mass hierarchy resolving power of the T2KK 
experiment, because the sensitivity grows as the difference in the magnitude of the earth 
matter effect between the oscillation probabilities observed at the near (SK) and a far 
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Figure 12: The same as Fig. HU but the 2.5° at SK. 

(Korea) detectors. 

As for the Tokai-to-Korea baselines, we find that the average matter density po grows 
from 2.85 g/cm 3 at L = 1000 km to 2.98 g/cm 3 at L = 1100 km, which are significantly 
higher than the PREM value of 2.70 g/cm 3 (L = 1000 km) and 2.78 g/cm 3 (L = 1100 
km); see Table [2J This is essentially because both the Conrad discontinuity between 
the upper and the lower crust, and the Moho- discontinuity between the crust and the 
mantle are significantly higher than PREM in the main region of the baseline below the 
Japan/East sea. This also contributes positively to the mass hierarchy measurement. 
The distribution is convex as in PREM but Repi/po is smaller in magnitude mainly 
becomes of the thin crust under the sea. One subtle point which we notice is that the 
L = 1000 km baseline almost touches the Moho-discontinuity at its bottom, see Fig. |5j 
and hence small error in the depth of the discontinuity can change the average density 
significantly. If the Moho-discontinuity at about 20 km below the sea level is shallower 
by 0.7 km, which is the present error of the seismological measurement, po grows from 
2.85 g/cm 3 to 2.90 g/cm 3 , by about 2%. On the other hand, this uncertainty does not 
affect the T2KK physics capability significantly because the major error in the earth 
matter density enters when the sound velocity data is converted to the matter density, 
for which we assign 6% error in our analysis; see Fig. [3j 

With the updated matter density distributions, we repeat the A% 2 analysis of Refs. 
[T?|[T8"]. which estimates the mass hierarchy resolving power of the T2KK experiment, by 
performing a thought experiment with a far detector of 100 kt fiducial volume in various 
location of Korea along the T2K beam line. In order to perform the full parameter scan 
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effectively, we introduce a new algorithm to compute the oscillation probabilities exactly 
for step-function-like matter density profile, where the earth matter effect is computed 
by diagonalizing real symmetric 3x3 matrix by decoupling it from the evolution due to 
the #23 mixing and the CP phase [34]. Analytic and semi-analytic approximations for 
the transition and survival probabilities are introduced to obtain physical interpretations 
of the numerical results. 

We confirm the previous findings of Refs. [TTJ [18] that the mass hierarchy resolving 
capability is maximized when the T2K beam is oriented downward to 3.0° OAB at SK 
and the far detector is placed in the south-east coast of Korean peninsula at L ~ 1000 km 
where ~ 0.5° OAB can be observed. In this report, we repeat the analysis by splitting 
the total beam time into half neutrino (v^) and half anti-neutrino (p^) beams, and find 
significant improvements in the mass hierarchy resolving power. Most remarkably, we 
find that the highest level of sensitivity can now be achieved for a far detector in locations 
where the T2K neutrino beam can be detected at an off-axis angle up to ~ 1°, because 
of the significant difference in the E v - dependence of the oscillation probability between 
Vfj, — > v e and v^ — > v e . This observation expands significantly the area in which a far 
detector can be most effective in determining the neutrino mass hierarchy pattern. 
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Appendix 

In the appendix, we present the second order perturbation formula of the time- 
evolution operator for the location (x) dependent Hamiltonian of eq. f lT5]) : 

H{x) = H + H X + 8H x {x) , (A.l) 

where each term has the matrix element 

2 2 

/ ,T7 I \ m2 2- m J T T TT* d ° X X 

(v p \ H 1 \u a ) = — — — U(3 2 U a2 + —dp e 8 ae , 



-1 oo 

(I/pi 8H x {x) \v a ) = -E 

h k=l L 



Re(afc) cos — - — J + Im(ofe) sin 





(A.2a) 


(A.2b) 


2nkx\ 
L J. 


5p e 8 ae .(A.2c) 



The first two terms, H and Hi, do not depend on the location x, whereas in the last 
term 8Hi(x) the x-dependent part of the matter density distribution along the baseline, 
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p(x) — po is expressed as a Fourier expansion. The second term in eq. ( 1A.2bj) gives the 
matter effect due to the average matter density p . 

In the leading order, we consider only Hq of eq. (1A.2a|) . The time evolution via Hq 
can readily be solved and the \v a ) — > \vp) transition matrix element is 



S (L) Pa = (u p \ e~ iH ° L \v a ) = U pi U* al + Up 2 U* a2 + Up 3 U* 3 e~ iA ™ 
= <W + Up 3 U* 3 (e- A - - l] 



(A.3) 



with A 13 as in eq. fl20l) . In the next order, we divide the corrections into two parts, 



Si{L)p a — S ± {L)p a + S x (L)p a , 
where the first term is linear in Hi, 



(A.4) 



7 (i) 



S[ '{L)p a = -i(vp\ 



f L dxe iHo{x - L) Hie~ iHoX J \v a ) 



-iA 12 [/g 2 [/* 2 - * 



2E 



5/3e5ae ~ Aa a + 2B ' a a 



+ (e- iAl3 -l)\B 0a +-£- (A? a - 2Bp a 



(A.5) 



and the second term is linear in 5 Hi. 



Si (L)j3a 



L e~ iAl3 - 1 
E A^ 

Re(a fc )A 2 



f L dxe iHo{x - L) 5Hi(x)e~ iHox ) \v a ) 



x 



E 
it=i 



Afg - (2nk) 



l:j (a 9R x Jm(o fc )27rfcA 13 



A? 3 - (2nk) 



T^Pa 



(A.6) 



Here the coefficients Ap a , Bp a , and Cp a are quadratic terms of the MNS matrix elements, 



Ap a = Sp e lI e3 U^ 3 + Up 3 U* 3 5 ae , 

Bp a = U/3 3 U^ 3 \U e3 \ , 

C/3a = OpJJ e3 \J* a3 — U/3 3 U* 3 5 ae . 



(A. 7a) 
(A.7b) 
(A.7c) 



Generally, the Bp a terms are much smaller than the other terms because of the constraint 

\U e3 \ 2 < 0.044. Since A^ = and |-B MM | ~ {U^U^l 2 <C 1, we can approximate S[ (L)^ 
rather accurately by — iA 12 |f/ M2 | 2 , which does not depend on the matter effect term a . 
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On the other hand, \A efl \ can be as large as \U v aU e 2\ in S\ (L) eM , and the matter effect 
can be significant. 

(2) 

As for the x-dependent perturbation S^ a , we note that the terms proportional to 
Re(afc) are in phase with the leading term of eq. ( 1A.3J) for (3a = e/jL, whereas those 
with Im(afc) are orthogonal. Therefore the Im(afc) terms contribute to the v^ — > v e 
oscillation probability only as Im(afc) 2 |C eAt | 2 ~ (a k ) 2 \U e3 \ 2 , which is strongly suppressed. 
Contributions to the v^ — > v^ survival probability are strongly suppressed again in this 



order due to the smallness of \A 



I'!' 



2B^\ and C^ 



0. 



Finally, the second order correction, S2, is evaluated only for x-independent perturba- 
tion, Hi, because the magnitude of all the non-zero Fourier coefficients are significantly 
smaller than the average, po; see Figs. HI |6]and Table [2j We find 



S 2 (L)p a = - (v p \ (£ dx £ dye lH ^- L ^H l e iH ^-^H 1 e 



i.H y 



\ U a) 



(A.* 



-i E I dx [-fU p2 U; 2 + || {5p e 5 ei + Up 3 U: 3 6 n (e^~^ - 1 



bi [X)j a 
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^U«,U 
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A12 apL 
2 2E 



/32^„ 2 
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2% 2 (e" A « - 1 
A13 A 2 3 



Ega 



OqL 

2E 



- {(l - \U e3 \ 2 ) (S $e 5 ae - A Pa ) + B Pa ) 



^13 



{<WU ^| 2 + (l - 3\U e3 \ 2 ) A 0a - 3 (l - 2\U e3 \ 2 ) B Pa ) 



+ (e^ Al3 - l) \^\U e3 \ 2 B Pa + -£- (\U e3 \ 2 Ap a + (l - 3\U e3 \ 2 ) B Pa ) 



*~£T {-^e5 ae \U e3 \ 2 - (l - 3\U e3 \ 2 ) A 0a + 3 (l - 2\U e3 \ 2 ) B pa 



,(A.9) 



where Dp a and Ep a are the MNS matrix element factors 

D/3a = Sp e U e 2Ua2 + U^U^ae , 

E Pa = u 02 u; 2 u e3 u* a3 + Up 3 u: 3 u e2 u* a2 . 



(A.lOa) 
(A.lOb) 



It should be noted that although A 2 2 is small, 0.01, around the first oscillation maximum 
1^13 1 ~ k, it is significant around the second oscillation maximum |A 13 | ~ 3tt where it 
grows to 0.1. The combination, A 12 (a L/2E) is about 0.02 for T2K and 0.06 for L ~ 1000 
km around |Ai3| ~ n. The matter effect term squared, (doL/2E) 2 , is 0.03 at Kamioka 
while it grows to ~ 0.3 or larger for the Tokai-to-Korea baselines. 
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